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Abstract 


Prediction  of  rigid  body  falling  through  water  column  with  a  high  speed  (such  as 
Mk-84  bomb)  needs  formulas  for  drag/lift  and  torque  coefficients,  which  depend  on 
various  physical  processes  such  as  supercavitation  and  bubbles.  A  diagnostic- 
photographic  method  is  developed  in  this  study  to  determine  the  drag/lift  and  torque 
coefficients  for  a  fast  moving  rigid-body  in  water  column.  The  diagnostic  part  is  to  derive 
the  relationships  (called  diagnostic  relationships)  between  (drag,  lift,  torque)  coefficients 
and  (position,  orientation)  of  the  rigid  body  from  the  three  momentum  equations  and 
three  moment  of  momentum  equations.  The  photographic  part  is  to  collect  data  of 
trajectory  and  orientation  of  a  fast  moving  rigid  body  using  multiple  high-speed  video 
cameras  (10,000  Hz).  Substitution  of  the  digital  photographic  data  into  the  diagnostic 
relationships  leads  to  semi-empirical  formulas  of  drag/lift  and  torque  coefficients,  which 
are  functions  of  the  Reynolds  number,  attack  angle,  and  rotation  rate.  This  method  was 
verified  by  1/12**'  Mk-84  bomb  strike  experiment  with  various  tail  configurations  (tail 
section  with  four  fins,  two  fins,  and  no  fin  and  no  tail  section)  conducted  at  the  SRI  test 
site.  The  cost  of  this  method  is  much  lower  than  the  traditional  method  using  the  wind 
tunnel.  Various  trajectory  patterns  are  found  for  different  tail  configurations. 


Keywords;  Body-flow  interaction,  drag  coefficient,  lift  coefficient,  torque  coefficient, 
fast  rigid  body  in  water  column,  dynamic-photographic  method 
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1.  Introduction 


Prediction  of  a  fast  moving  rigid  body  in  water  column  has  wide  scientific 
significance  and  technical  application.  The  dynamics  of  a  rigid  body  allows  one  to  set  up 
six  nonlinear  equations  for  the  most  general  motion:  three  momentum  equations  and  three 
moment-of-momentum  equations.  The  scientific  studies  of  the  hydrodynamic 
characteristics  of  a  rigid  body  in  water  column  involve  nonlinear  dynamics,  body  and 
multi-phase  fluid  interaction,  supercavitation,  and  instability  theory. 

The  technical  application  of  the  hydrodynamics  of  a  rigid  body  with  high  speed 
into  fluid  includes  aeronautics,  navigation,  and  civil  engineering.  Recently,  the  scientific 
problem  about  the  movement  of  a  rigid  body  in  water  column  drew  attention  to  the  naval 
research.  This  is  due  to  the  threat  of  mine  and  maritime  improvised  explosive  device 
(lED).  Prediction  of  a  fast  falling  rigid  body  in  the  water  column  contributes  to  the  bomb 
breaching  for  mine/maritime  lED  clearance  in  surf  and  very  shallow  water  zones  depth 
shallower  than  12.2  m  (i.e.,  40  ft),  according  to  El.S.  Navy’s  standards. 

To  predict  rigid  body  maneuvering  in  water  column,  accurate  calculation  of 
hydrodynamic  forces  and  torques  is  important  but  difficult.  Difficulty  is  caused  by 
unknown  drag,  lift,  and  torque  coefficients.  Thus,  determination  of  these  coefficients  is 
crucial  for  the  prediction'^ Different  from  traditional  methods  used  in  aerodynamics 
such  as  using  the  wind  tunnel,  we  present  a  new  efficient  and  low  cost  method  to 
determine  the  drag,  lift,  torque  coefficients.  This  method  consists  of  dynamical  and 
photographical  parts.  The  dynamical  part  is  to  establish  diagnostic  relationships  between 
hydrodynamic  coefficients  (e.g.,  drag,  lift,  torque  coefficients)  and  kinematics  of  the  rigid 
body  (e.g.,  position,  orientation).  The  photographical  part  is  to  conduct  experiments  with 
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shooting  the  rigid  body  into  the  hydrographieal  tank  with  high  speed  and  reeording  its 
underwater  position  and  orientation  with  high-speed  video  (HSV)  eameras  at  10,000 
films  per  second  (fps).  Semi-empirical  formulas  can  be  established  for  the  drag,  lift, 
torque  coefficients  by  substituting  the  digital  data  of  rigid  body’s  location  and  orientation 
from  the  HSV  cameras  into  the  diagnostic  relationships. 

The  rest  of  the  paper  is  outlined  as  follows.  Sections  2  and  describes 
hydrodynamic  forces/torques  and  basic  dynamics.  Section  4  presents  the  diagnostic 
relationships  for  the  drag/lift  and  torque  coefficients,  which  were  derived  from  the 
momentum  and  moment  of  momentum  equations.  Section  5  depicts  the  photographic 
method.  Sections  6  and  7  show  the  semi-empirical  formulas  for  the  drag/lift  and  torque 
coefficients  and  the  verification.  Section  8  describes  the  tail  section  damage  effects. 
Section  9  presents  the  conclusions. 

2.  Hydrodynamic  Forces  and  Torques 

2,1.  Drag/Lift  Forces 

Let  the  earth-fixed  coordinate  system  is  use  with  (i,  j)  the  unit  vectors  in  the 
horizontal  plane  and  k  the  unit  vector  in  the  vertical  direction.  Consider  an  axially 
symmetric  rigid  body  such  as  a  bomb  falling  through  water  column.  The  two  end-points 
of  the  body  (i.e.,  head  and  tail  points)  are  represented  by  r/,(t)  and  r^O-  The  difference 
between  the  two  vectors  in  nondimensional  form 

(1) 

is  the  unit  vector  representing  the  body’s  main  axis  direction  (Fig.  1).  The  centers  of  mass 
{om)  and  volume  (ov)  are  located  on  the  main  axis  with  cr  the  distance  between  Ov  and 
Om,  which  has  a  positive  (negative)  value  when  the  direction  from  Ov  to  Om  is  the  same 
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(opposite)  as  the  unit  veetor  e  (Fig.  2).  The  loeation  (or  ealled  translation)  of  the  body  is 
represented  by  the  position  of  Om, 

r(0  =  x\  +yj  +  zk.  (2) 

The  translation  veloeity  is  given  by 

dv 

-^  =  n,  u  =  Ue^,  (3) 

dt 


where  {U,  e„)  are  the  speed  and  unite  veetor  of  the  rigid-body  veloeity.  Let  jd  be  the 
angle  of  the  body’s  main  axis  with  the  horizontal  plane,  ybe  the  angle  of  the  body’s 
veloeity  u  with  the  horizontal  plane,  and  a  be  the  attaek  angle  between  the  direetion  of 
the  main  body  axis  (e)  and  the  direetion  of  the  body  veloeity  (Cu)  (Fig.  2), 

a  =  eos  '  (e„*e),  =  sin  '  (e*k),  y  =  sin  '  (e„*k) .  (4) 


Usually,  the  unit  veetor  e„  is  represented  by  (Fig.  3) 

e„  =  eos  y  eos  +  eos  y  sin  +  sin  yk  . 

where  y/  is  the  azimuth  angle.  Differentiation  of  (5)  with  respeet  to  t  gives 


dQ 


dy  ,  dy/ 


-  =  — ^e„  +^-eosye 
dt  dt  dt 


V 

U  ’ 


where 


= -sin^i  +  eos^^j,  =-sinycos^i-sinysin^j-i-eosyk. 


are  unit  vectors.  It  is  clear  that 


(5) 


(6) 


(V) 


e*"  1  e  ,  e"  1  e  ,  e"  1  e*" . 

U  ll  ^  u  u  ^  u  u 


(8) 


Let  Vw  be  the  water  velocity.  Water-to-body  relative  velocity  V  (called  the 
relative  velocity)  is  represented  by 

V^V„-u«-u  =  -Ue„.  (9) 
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Here,  the  water  veloeity  is  assumed  mueh  smaller  than  the  rigid-body  velocity.  A  third 


basic  unit  vector  ( e* )  can  be  defined  perpendicular  to  both  e  and  e^, 


m  I  I 

e.,  xe 


(10) 


The  drag  force  (F^)  is  in  the  opposite  direction  of  the  rigid-body  velocity.  The  lift  force 
(F/)  is  in  the  plane  constructed  by  the  two  vectors  (e,  e^)  (i.e.,  perpendicular  to  e* )  and 


perpendicular  to  e„,  and  therefore  the  lift  force  is  in  the  direction  of  e*  x  e„ .  Both  drag 


and  lift  forces,  exerting  on  the  center  of  volume,  are  represented  by 

“  ~fd^u  ’ 


(11) 


where  fd  and  fi  are  the  magnitudes  of  the  forces.  The  magnitudes  {fd,  fi)  are  represented 
by  the  drag  law^^^. 


(12) 


where  p  is  water  density;  is  the  under  water  projection  area;  (C*  C/)  are  the  drag 
and  lift  coefficients,  which  are  determined  by  the  experiments. 

Bomb  is  a  fast  moving  rigid-body  usually  with  four  fins.  Two  fins  in  the  same 
plane  are  called  the  pair  of  fins.  For  simplicity,  these  fins  are  treated  approximately  as 
the  NACA0015  airfoils  (Fig.  4)  with  known  drag  and  lift  coefficients.  For  example,  they 
are  listed  at  .  Using  these 

coefficients,  the  combined  drag  and  lift  forces  on  a  pair  of  fins  (  F/  )  can  be  calculated.  If 
the  bomb  has  n  pairs  of  fins,  the  total  drag  and  lift  forces  on  the  fins  are  represented  by 
«F/. 


2,2,  Translation  Torque 
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As  the  drag/lift  forces  are  defined  on  the  center  of  volume,  the  hydrodynamic 
torque  (relative  to  center  of  mass)  M/,  is  caused  by  (a)  drag  and  lift  forces  exerted  on  Ov 
to  Om,  (b)  moment  of  drag/lift  forces  exerted  on  Ov  (Mtr)  and  (c)  anti-rotation  of  the  rigid- 
body  (Mrot), 

M,  =-o-ex(F^+F;)-no-^exF/'  +  .  (13) 

The  torque  Mtr  is  perpendicular  to  both  e„  (the  direction  of  V)  and  e  (the  body 
orientation),  and  therefore  it  is  in  the  same  direction  of  the  unit  vector  e*  , 


Mtr  =  MX, 


(14) 


where  Mtr  is  the  magnitude  of  the  anti-translation  torque,  and  be  calculated  by  the  drag 
law'-'^^ 

U„=\c„pA,LJj\  (15) 


where  Cm  is  the  anti-translation  torque  coefficient. 

2,3.  Anti-Rotation  Torque 

The  anti-rotation  torque  acts  as  the  rigid-body  rotates.  Let  be  the  rigid-body’s 
angular  velocity  vector,  which  is  decomposed  into  two  parts  with  one  along  the  unit 
vector  e  (bank  angle)  and  the  other  H  (azimuthal  and  elevation  angles)  perpendicular  to 
e(Fig.  5), 

n*  =  Q^e  +  n.  (16) 


Let  e  be  the  unit  vector  in  the  direction  of  Cl , 

(O  ‘ 

n  =  Qe^,  Q  =  |n|.  (17) 

The  unit  vector  is  perpendicular  to  e. 
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(18) 


e,  *6  =  0 . 


Time  rate  of  ehange  of  the  unit  veetor  e  (main  axis  direetion)  is  given  by 

—  =  n*xe  =  nxe.  (19) 

dt 

Veetor  produet  between  (19)  and  the  unit  veetor  e  gives 

d^ 

ex  — =  n  =  Qe^.  (20) 

dt 

Differentiation  of  (16)  with  respeet  to  time  t  and  use  of  (20)  lead  to^^^ 

- =  — ^e  +  Q  (nxe)H - .  (21) 

dt  dt  ’  dt 


The  anti-rotation  torque  (Mm^)  is  against  the  rotation  of  the  rigid  bodyll* .  It  ean 
be  deeomposed  into  two  parts, 

M,,,=M,+M,,  (22) 


where  the  torque  (resistant  to  self  spinning,  Q^e  )  parallels  the  main  axis  of  the 
body  (i.e.,  the  unit  veetor  e), 

M^=-M^e,  (23) 


and  the  torque  is  perpendieular  to  the  unit  veetor  e. 


e^Te,  (24) 

where  Ms  and  Me  are  the  eorresponding  sealar  parts.  The  drag  law  shows  that'-’^ 

u,  =  \c,pAXp.\n.,  (25) 

M,  =  i  CF(ClL,  /  F,  )pA.LJ^  .  (26) 
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where  Cg  is  the  torque  eoeffieient  for  self-spinning;  C  is  the  drag  eoeffieient  due  to 
eross-body  flow.  The  funetion  F{ju)  is  obtained  from  surface  integration  of  torque  due  to 
cross-body  hydrodynamic  force  (perpendicular  to  the  body).  For  a  cylindrical  body,  the 
coefficient  C  and  function  F{ju)  are  given  by^^^ 


C  = 


1.9276  +  8/Re, 
1.261+16/Re, 
0.855+89/Re, 
0.84+0.00003Re, 
1.2-4/^ 
0.835  -  0.35/^, 
0.7-0.08/^, 


1.875 -0.0000045  Re, 
1/(641550/Re+  1.5), 


if  Re  <  12 
ifl2<Re  <180 
if  180  <  Re  <2000 
if2000<  Re  <12000 
if  12000  <  Re  <150000,  ^>10 
if  12000  <  Re  <150000,  2  <^<10 
if  12000  <  Re  <150000,  S<2 
if  150000  <  Re  <350000 
if  Re  >  350000. 


F{f)^ 


J_ 

6ju 


1  2 
- F 

v4  y 


4  2  1 

H - JU  H - :r 

3  2//" 


M  4^ 

- F 

vl6  y 


for  //  >  1  /  2 
for  //  <  1  /  2 


(27) 


(28) 


Here,  S  is  the  aspect  ratio  of  the  rigid  body;  Re  =  UD  /v ,  is  the  Reynolds  number  with 
D  the  effective  diameter  of  rigid  body;  is  the  projection  of  the  water-to-body  relative 
velocity  on  the  vector  x  e , 


n=V.e,=-f/e„.(e^xe),  (29) 

where  Eq.(9)  is  used.  For  a  supercavitation  area,  a  correction  factor  may  be  multiplied  to 
the  coefficient  C. 

3.  Dynamics 

3.1.  Momentum  Equation 
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Differentiation  of  (3)  with  respect  to  time  t  gives  the  acceleration  of  the  rigid 


body, 


dn  dU  ,rde„ 
—  =  — e„  +U — 
dt  dt  dt 


(30) 


The  momentum  equation  in  the  earth-fixed  coordinate  system  is  given  by  (Fig.  3), 


m 


dU 

dt 


t.+U 


dt 


=  F. 


(31) 


where  F  is  the  resultant  force  consisting  of, 


F  =  F^+F,+F,+F,+nF/ 


(32) 


Here, 

Fg=-mgk,  F^=/?ngk,  (33) 

are  the  gravity  and  buoyancy  force  and  Tt  is  the  volume  of  the  rigid  body.  For  bomb 
without  fin  {n  =  0),  the  resultant  force  is  represented  by 

F  =  F^+F,  +  F,+F,.  (34) 


Inner  products  between  Eq.(3 1)  and  the  unit  vectors  ( e„,  e^,  )  for  n  =  0  lead  to 


m— =  (/?n-m)gk.e„-/^,  (35) 

dt 

mU^  =  [{pIl-m)g^  + (36) 
mf/cosy^^  =  [(yC>n-m)gk-i-/;e,]»eJ^.  (37) 


Here,  we  use  (8)  and  the  condition  T  (i.e.,  the  lift  force  perpendicular  to  the  drag 
force).  The  vector  [(pH  -  m)gk  +  /^e,]  in  (36)  and  (37)  can  be  represented  by 


[(pH  -  m)gk  +  /,ej  =  {m^  +  /Je„  +mU^tl+mU  cosy  ,  (38) 

dt  dt  dt 
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which  shows  that  this  vector  has  the  component  of  (mdU/dt  +  fd)  in  the  direetion  of  e„, 
mUdy ! dt  in  the  direction  of  ,  and  mU cosydy/ ! dt  in  the  direction  of  .  Inner 
product  of  (38)  by  the  unit  vector  e/  leads  to 


[(-m  +  /?n)gk  +  /,ej»e; 


r  TT  dy/  ^ 

mu  — e  +mU  cosy — e 
dt  "  dt  “ 


•e 


I  ’ 


(39) 


where  the  condition  (e^  _L  e, )  is  used  [see  (11)]. 

3.2.  Moment  of  Momentum  Equation 

The  moment  of  momentum  equation  (relative  to  center  of  mass)  is  given  by 

J.^  =  M.+M.,  (40) 

dt 

where 


M^  = -ere  X  (pngk) ,  (41) 

is  the  torque  eaused  by  buoyancy  force  (called  the  buoyancy  torque).  J  is  gyration  tensor. 
In  the  body-fixed  eoordinate  system  J  is  a  diagonal  matrix  [7] 


J  = 


0 

0 


0  0 
J,  0 
0 


(42) 


with  J\  ,  J2,  and  J3  are  the  moments  of  inertial.  For  axially  symmetrie  rigid-body  sueh  as 
bomb,  J2=J2.  Substitution  of  (13)  into  (40)  and  use  of  (1 1)  and  (22)  lead  to 

j.^  =  M,+M,  (43) 

dt 


where 


Ms-cr/7ngexk-cr(/^exe„+  y]exe,)-u<7^//exef  +M,^-tM^.  (44) 
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Since  is  the  anti-self-spinning  torque  whieh  parallels  the  unit  veetor  e,  and  M  is  the 
torque  perpendieular  to  the  unit  veetor  e,  the  moment  of  momentum  equation  (43)  ean  be 
deeomposed  into  two  eomponents  with  one  parallel  to 


_ 

dt 


(45) 


and  the  other  perpendieular  to  e, 


da  M  ,  . 

-=--n.(nxe), 


(46) 


where  (24)  is  used.  For  a  rigid-body  with  very  slow  or  no  self  spinning  (sueh  as  bomb), 
==  0,  Eq.(46)  beeomes 


.  da  - 

J, - =  M  . 

'  dt 


(47) 


4.  Diagnostic  Relationships 


The  momentum  equations  (35)  and  (39)  ean  be  rewritten  by 


_2[(^n  -  m)  gk*e„  -  mdU  /  Jt] 


(48) 


mU{^eld/ / dt  +  el  cos /di// / dtye,  -[pTl-m 


(49) 


Inner  produets  of  (47)  by  the  unit  veetors  e*  (direetion  of  M^^)  for  without  fm  {n  =  0)  and 
(direetion  of  M^g)  lead  to 


da 


‘^2^*el  +  ^png(exk).e^ 


dt 


■  +  CF(Q4/FJ^e^.el 


(50) 


+^[Q  (e  X  >el  +  C,  (e  x  e,).e^] , 
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Eqs.(48)-(50)  are  the  diagnostic  relationships  for  determining  the  coefficients  (C*  C/, 
Cm).  In  the  righthand  sides  of  (48)-(50),  the  unit  vectors  (e,  Cv,  e^)  are  essential,  and  the 

other  unit  vectors  (e/,  e*  ,  e^,  )  are  derived. 

Accurate  prediction  of  a  high-speed  rigid  body’s  location  and  orientation  in  the 
water  column  needs  realistic  values  of  the  drag/lift  and  torque  coefficients  (C*  C/,  Cm). 
Among  these  coefficients,  (Q,  C/,  Cm)  depend  on  the  attack  angle  (a). 

Effects  of  special  phenomena  such  as  surface  impact,  bubbles,  and  supercavition 
on  the  movement  of  rigid  body  are  represented  in  the  change  of  these  coefficients.  Thus, 
if  the  time  evolutions  of  unit  vectors  (e,  tu,  e^)  and  variables  {x,y,z,U,y/,x)  are 

measured,  time  evolutions  of  the  drag/lift  and  torque  coefficients  (Cj,  C/,  Cm)  can  be 
obtained  using  the  diagnostic  relationships  (48)-(50).  The  rotation  rate  fl  is  calculated 
from  the  time  series  of  ( 

With  large-amount  of  derived  (C*  C/,  Cm)  data,  instantaneous  relationships  (semi- 
empirical  formulas)  can  be  established  statistically  between  (Q,  C/,  Cm)  and  basic 
parameters  such  as  the  attack  angle  a  ,  Reynolds  number  (Re),  and  the  rotation  rate  fl . 
A  traditional  method  for  this  purpose  is  to  conduct  experiments  in  the  wind  tunnel.  Else 
of  wind  tunnel  may  be  feasible  for  determining  the  drag/lift  and  torque  coefficients  of  a 
rigid  body  in  the  air,  but  not  realistic  in  the  water  especially  the  rigid  body  with  high 
speed  such  as  bomb. 

5.  Photographic  Method 

5.1.  General  Description 

As  alternative  to  the  wind  tunnel  method,  the  photographic  method  is  used  to 
measure  time  evolutions  of  rigid-body’s  head  and  tail  points  rk(t)  and  r,(0  using  HSV 
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cameras 


From  the  [r/,(^),  r,(0]  data,  the  time  series  of  the  unit  veetors  (e,  Cj,,  e^), 
variables  {x,y,z,U,y/,x),  and  attack  angle  («)  can  be  calculated  using  (3),  (4),  (5), 
and  (23).  With  the  diagnostic  relationships  (48)-(50),  time  evolutions  of  the  drag/lift 
and  torque  eoefficients  (Cj,  C/,  Cm)  can  be  obtained. 

5.2,  Model  Design  and  Construction 

Details  of  the  full-seale  Mk84  with  a  guidance  tail  seetion  were  obtained  from 
Boeing  Corporation  with  the  mass  inertia  properties  for  the  eomplete  Mk84  warhead  and 
tail  section  as  well  as  for  the  modifications  to  the  tail  section  including  removal  of  two 
fins,  four  fins,  and  the  complete  tail  seetion.  We  developed  a  closed  form  analytical 
model  to  determine  the  outer  casing  material  and  thickness  and  type  of  internal 
components  and  their  loeation  for  the  1/12-scale  model  to  mateh  the  scaled  total  mass, 
eenter  of  gravity  and  radius  of  gyration  values  for  the  Mk84  bomb  with  the  different  tail 
eonfigurations.  Here,  the  overall  outer  shape  of  the  scaled  Mk84  bomb  represented  an 
exact  dimensional  replica  of  the  full-scale  system.  Fig.  6  shows  the  overall  design  details 
of  each  model  configuration.  The  outer  casing  was  made  from  7075-T6  aluminum.  For 
the  models  with  a  tail  seetion,  the  easing  was  fabricated  in  two  pieces  that  were  screwed 
together  near  the  center  point.  To  obtain  the  eorreet  mass  properties,  a  copper  plug  was 
inserted  inside  the  model  and  then  the  remaining  internal  eavity  was  filled  with  a  low 
density  Epoxy.  For  all  of  the  models  with  a  tail  section,  only  small  differences  in  the 
copper  plug  size  was  required  to  match  the  mass  properties.  For  these  models  the  Epoxy 
had  a  density  of  0.546  kg  m'  .  Eor  the  model  with  no  tail  seetion,  two  eopper  plugs  were 
required  and  the  Epoxy  density  had  to  be  inereased  to  1.168  kg  m'  by  adding  Tungsten 
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powder.  Four  types  of  model  Mk-84  bombs  were  constructed;  Type-I  (tail  with  4  fins), 
Type-II  (tail  with  two  fins),  Type-Ill  (tail  with  no  fins),  and  Type-IV  (no  tail). 

5,3,  Experiment  Procedure 

Two  experiments  were  conducted  to  demonstrate  the  feasibility  of  this  method 
with  1/12*''  scale  model  of  the  general  purpose  bomb  (Mk-84)  as  the  fast  moving  rigid 
body.  The  first  experiment  was  conducted  at  the  Monterey  Bay  Aquarium  Research 
Institute  (MBARI)  Unmanned  Underwater  Vehicle  Test  Tank  (9.14  m  x  13.72  m  x  9.14 
m)  filled  with  standard  sea  water^'^^'^'^^.  The  second  experiment  was  conducted  at  a  6  m 
deep  by  9  m  diameter  pool,  located  at  the  SRI’s  Corral  Hollow  Experiment  Site'-'"'l  A 
pneumatic  launcher  (gas  gun)  was  used  to  shoot  the  1/12"'  model  Mk-84  bomb  into  the 
water  tank  in  the  first  (second)  experiment  with  velocities  up  to  100  m  s"'  (454  m  s"')  (Fig. 
7).  Since  the  implementation  is  very  similar  between  the  two  experiments  and  the  second 
experiment  has  much  larger  speed,  the  second  experiment  is  presented  here  to  illustrate 
the  photographic  method. 

Models  of  Mk-84  bombs  with  and  without  tail  section  are  taken  as  examples  to 
illustrate  the  methodology  for  determination  of  the  drag/lift  and  torque  coefficients  (C* 
Cl,  Cm),  and  in  turn  the  prediction  of  location  and  orientation  of  a  fast-moving  rigid-body 
through  the  water  column.  The  primary  objective  is  to  determine  the  Mk84  trajectory 
through  the  very  shallow  water  zone  to  provide  an  estimate  of  the  maximum  bomb-to- 
target  standoff  and  required  fuse  delay  time  for  optimum  target  lethality.  Because  it  is 
possible  that  a  portion,  or  all,  of  the  guidance  tail  section  may  become  separated  from  the 
warhead  during  water  entry,  it  is  necessary  to  determine  the  Mk84  trajectory  for  a  variety 
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of  different  tail  eonfigurations  ranging  from  a  warhead  with  a  completely  intact  tail 
section  and  four  fins  to  a  warhead  with  the  tail  section  completely. 

Using  the  Hopkinson  scaling  laws,  1/12-scale  Mk84  bomb  models  were  designed 
and  constructed  in  SRI  that  matched  the  overall  casing  shape  and  mass  inertial  properties 
of  the  full-scale  Mk84  prototype.  To  model  the  different  possible  damaged  tail 
configurations,  we  fabricated  models  that  consisted  of  the  warhead  section  with  a 
complete  tail  section  and  four  fins,  a  complete  tail  section  and  two  fins,  a  complete  tail 
section  and  no  fins,  and  with  the  tail  section  removed. 

The  gas  gun  (0.10  m  diameter  and  1.52  m  long)  barrel  was  evacuated  before 
launching  the  scale  model  to  prevent  an  air  blast  from  disturbing  the  water  surface  prior 
to  the  model  impacting  the  water  surface.  At  the  end  of  the  gas  gun  there  was  a  massive 
steel  ring  to  strip  the  sabot  from  the  scale  model.  At  high  velocities  there  is  some 
deviation  from  the  theoretical  calibration  curve,  which  may  be  attributed  to  gas  blow  by 
around  the  sabot  or  friction.  For  the  maximum  gun  operating  pressure  of  2,500  psi,  we 
were  able  to  achieve  a  nominal  water-entry  velocity  of  about  305  m  s"'. 

Two  orthogonal  periscope  housings  were  positioned  in  the  pool  to  allow 
simultaneous  above-water  and  below-water  visualization  of  the  model  trajectory.  The 
housings  supported  Phantom-7  HSV  cameras,  which  were  run  at  10,000  fps.  Five  high- 
intensity,  short  duration  (30  ms)  flash  bulbs  were  used  to  front-light  the  scale  model  as  it 
entered  the  water  and  traveled  under  water.  The  HSV  cameras  and  flash  bulbs  were 
triggered  at  the  time  the  sabot  was  released  within  the  gun. 

5.4,  Data  Collection 
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For  the  complete  Mk84  bomb  system,  including  the  warhead  with  tail  section  and 
4  fins  (Type-I),  tail  section  and  2  fins  (Type-II),  tail  section  and  no  fin  (Type-Ill),  and  no 
tail  section  (Type-IV)  were  launched  at  different  nominal  water-entry  velocity  regimes 
from  119.48  m  s"'  to  308.83  m  s'\  Table  1  summarizes  the  overall  experimental  matrix 
and  water-entry  conditions.  Typically,  the  water-impact  angle  of  entry  was  between  88° 
and  90°.  In  Launches  10,  II,  and  12  the  sabot  failed  to  fully  support  the  scale  model 
within  the  gun  during  the  launch  phase,  resulting  in  the  scale  model  impacting  the  sabot 
stripper  plate  before  impacting  the  water.  Sequences  of  images  from  the  two  orthogonal 
HSV  cameras  were  generated  for  each  launch  (Fig.  8). 

Therefore,  16  time  series  of  the  unit  vectors  (e,  tu,  e^),  variables 
(x,y,z,17,^,;j',n  ),  and  attack  angle  («)  were  obtained  for  the  four  types  of  the 
model  Mk-84  bombs.  All  the  experimental  data  have  been  converted  to  full-scale  values. 
We  divided  the  16  time  series  into  two  groups:  (a)  Launches  13,  14,  15  (called  the 
working  data)  for  determining  semi-empirical  formulas  for  the  drag/lift  and  torque 
coefficients  (Cj,  C/,  Cm),  and  (b)  rest  of  the  data  for  evaluating  the  semi-empirical 
formulas  (called  evaluation  data). 

6.  Semi-Empirical  Formulas 


Statistical  analysis  was  conducted  on  the  working  dataset  (Type-IV:  Launch- 1 3,  - 
14,  and  -15)  at  each  time  step  (remember  the  data  rate  is  10,000  Hz)  between  (C*  C/,  Cm) 
and  (Re,  a  ,  fl).  Following  semi-empirical  formulas  have  been  established. 


C  =0.02  +  0.35e 


f  Re  ^ 


0.2 


vRe*y 


-I-  0.008Qsm6* , 


(51)  (53) 
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Here,  Re*  =  1.8  x  10^,  is  the  eritieal  Reynolds  number,  and 

1 

6  =  -[n -\n -2a)  ,  (54)  (56) 

(2aV  (2a  T 

e,=7t  —  ,  6»2=2;r  — -1  .  (55)  (57) 

\  7t  )  \  7t  ) 

The  semi-empirieal  formulas  (51)-(53)  show  that  the  drag/lift  eoeffieients  (Q,  C/)  depend 
more  on  (Re,  a)  and  less  on  the  rotation  rate  Q .  For  the  same  attaek  angle  («),  Cd 
inereases  with  Re.  For  the  same  Re,  Cd  inereases  with  a  monotonieally  from  0“  to  90° 
and  reduees  monotonieally  with  a  from  90°  to  180°  with  a  maximum  value  for  a  =  90° 
(Fig.  9).  The  dependenee  of  lift  eoeffieient  C/  on  (Re,  « )  is  a  little  eomplieated  than  Q, 
espeeially  for  the  attaek  angle  larger  than  90°  (Fig.  10).  The  torque  eoeffieient  Cm 
depends  only  on  (Re,  a )  (Fig.  1 1). 

The  momentum  equations  (35)-(37)  and  the  moment  of  momentum  equation  (47) 
were  integrated  using  the  same  parameters  (sueh  as  the  density  ratio,  length,  radius,  the 
eenter  of  volume,  and  the  eenter  of  mass)  and  the  drop  initial  eonditions  (speed  and 
orientation)  as  in  Type-111  (tail  without  fin)  and  Type-lV  (no  tail)  (see  Table  1)  after 
using  (51)-(53)  for  (Q,  C/,  Cm).  The  validity  of  the  semi-empirieal  formulas  (51)-(53) 
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are  verified  through  eomparison  between  ealeulated  and  observed  bomb  trajeetories, 
orientations,  and  veloeities. 

7.  Results  and  Verification 

7.1.  Type-I 

This  type  is  for  the  eomplete  Mk84  bomb  system,  ineluding  the  warhead  with  tail 
seetion  and  four  fins.  Seven  launehes  (Launeh-1  to  Launoh-7)  were  eondueted  at  different 
nominal  water-entry  veloeity  regimes  (119  ms"'  to  302  ms'^).  Fig.  8  shows  a  sequenee  of 
images  from  the  two  orthogonal  HSV  eameras  for  a  nominal  veloeity  of  295  ms"' 
(Launoh-3).  The  eavitated  eolumn  generated  by  the  Mk84  bomb  motion  eonsists  of  a 
tapered  eone  that  has  a  maximum  full  seale  diameter  at  the  end  of  the  tail  seetion  of  about 
0.9  m.  This  is  about  a  faetor  of  two  larger  than  the  maximum  bomb  diameter  of  0.45  m. 
This  eavity  shape  was  about  the  same  for  all  of  the  initial  water-entry  veloeities  between 
119  ms'^  and  302  ms'\ 

Fig.  12  shows  the  eomparison  of  bomb’s  translation  and  orientation  for  Type-I 
(Launeh-1  to  Launoh-7)  between  the  ealeulated  and  observed  data.  Both  ealeulated  and 
observed  trajectories  show  similar  patterns.  For  the  low  velocity  regime  of  about  125  ms"' 
(Launach-1  and  Launch-7),  at  a  full  scale  depth  of  12  m,  the  horizontal  position  ranged 
between  0.1  m  (Launch-7,  Fig.  12g)  and  0.67  m  (Launch-1,  Fig.  12a).  The  bomb 
trajectories  are  quite  stable  without  oscillation  and  tumbling  no  matter  the  water  entry 
velocity  is  high  or  low. 

7.2.  Type-II 

This  type  is  for  the  modified  Mk84  bomb  system  including  the  warhead  with  a  tail 
section  and  two  fins.  Three  launches  (Launch- 10,  -11,  -19)  were  conducted  at  an  average 
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water-entry  veloeity  of  about  294  ms''.  Fig.  13  shows  a  sequenee  of  images  (Launeh-1 1) 
from  the  two  orthogonal  HSV  eameras  with  a  water  entry-veloeity  of  290  ms"'.  The 
eavitated  eolumn  generated  by  the  bomb  motion  eonsists  of  a  tapered  eone  that  has  a 
maximum  full-seale  diameter  at  the  end  of  the  tail  seetion  of  about  0.9  m.  This  is  about  a 
faetor  of  two  larger  than  the  maximum  bomb  diameter  of  0.45  m.  Thus,  the  initial  eavity 
shape  was  about  the  same  as  for  the  model  with  a  tail  seetion  and  four  fins  as  deseribed 
above. 

Fig.  14  shows  the  eomparison  of  bomb’s  translation  and  orientation  for  Type-11 
(Launeh-10,  -11,  -19)  between  the  ealeulated  and  observed  data.  Both  ealeulated  and 
observed  trajeetories  show  similar  patterns.  At  full-seale  depth  of  12.2  m  (i.e.,  40  ft),  the 
horizontal  position  ranged  between  0.53  m  (Launeh-1 1)  and  2.1  m  (Launeh-10).  These 
values  are  about  a  faetor  of  two  larger  than  the  values  measured  for  Mk84  bomb 
eonliguration  with  a  tail  seetion  and  four  fins.  Also,  there  seems  to  be  no  eorrelation 
between  trajeetory  path  and  initial  impaet  angle.  The  removal  of  two  fins  eauses  the 
bomb  to  eventually  makelSO”  turn  and  travel  toward  the  surfaee.  The  12.2  m  depth  is 
reaehed  at  about  45  ms  after  water  entry. 

7.3.  Type-III 

This  type  is  for  the  modified  Mk84  bomb  system  ineluding  the  warhead  with  a  tail 
seetion  and  no  fin.  Three  launehes  (Launeh-1 6,  -17,  -18)  were  eondueted  at  an  average 
water-entry  veloeity  of  about  298  ms''.  Fig.  15  shows  a  sequenee  of  images  (Launeh-1 7) 
from  the  two  orthogonal  HSV  eameras  with  a  water  entry-veloeity  of  297  ms''.  The 
eavitated  eolumn  generated  by  the  bomb  motion  eonsists  of  a  tapered  eone  that  has  a 
maximum  full-seale  diameter  at  the  end  of  the  tail  seetion  of  about  0.9  m.  This  is  about  a 
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factor  of  two  larger  than  the  maximum  bomb  diameter  of  0.45  m.  Thus,  the  initial  eavity 
shape  was  about  the  same  as  for  the  model  with  a  tail  seetion  and  four  fins  and  for  the 
model  with  two  fins  as  deseribed  above. 

Fig.  16  shows  the  eomparison  of  bomb’s  translation  and  orientation  for  Type-Ill 
(Launeh-16,  -17,  -18)  between  the  ealeulated  and  observed  data.  Both  ealeulated  and 
observed  trajeetories  show  similar  patterns.  At  full-seale  depth  of  12.2  m  (i.e.,  40  ft),  the 
horizontal  position  ranges  between  2.1  m  (Launeh-18)  and  3.5  m  (Launeh-17).  These 
values  are  about  a  faetor  of  9.5  and  3.5,  respeetively,  larger  than  the  values  measured  for 
a  Mk84  warhead  with  a  tail  seetion  and  four  fins.  Also,  there  seems  to  be  no  eorrelation 
between  trajeetory  path  and  angle  of  impaet.  As  shown  in  Figs.  15  and  16,  the  removal  of 
four  fins  eauses  the  bomb  eventually  to  make  180'’  turn  and  travel  toward  the  surfaee  in  a 
manner  similar  to  the  trajeetory  for  a  bomb  with  a  tail  seetion  and  two  fins.  The  12.2  m 
depth  is  reaehed  at  about  47  ms  after  water  entry.  In  general,  the  model  with  a  tail  seetion 
and  no  fins  shows  a  deerease  in  overall  trajeetory  stability  eompared  with  a  model  with  a 
tail  seetion  and  two  fins.  This  is  evident  primarily  through  larger  horizontal  positions. 

7.4,  Type-IV 

This  type  is  for  the  modified  Mk84  bomb  system  ineluding  the  warhead  and  no 
tail  seetion.  Three  launehes  (Launeh-13,  -14,  -15)  were  eondueted  at  an  average  water- 
entry  veloeity  of  about  299  ms"'.  Fig.  17  shows  a  sequenee  of  images  (Launeh-13)  from 
the  two  orthogonal  HSV  eameras  with  a  water  entry-veloeity  of  297  ms'*.  The  eavitated 
eolumn  generated  by  the  bomb  motion  eonsists  of  a  tapered  eone  that  has  a  maximum 
full-seale  diameter  at  the  end  of  the  tail  seetion  of  about  0.9  m.  This  is  about  a  faetor  of 
two  larger  than  the  maximum  bomb  diameter  of  0.45  m.  Thus,  the  initial  eavity  shape 
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was  about  the  same  as  for  the  model  with  a  tail  section  (four,  two,  and  no  fins)  as 
described  above. 

Fig.  18  shows  the  comparison  of  bomb’s  translation  and  orientation  for  Type-IV 
(Launch- 13,  -14,  -15)  between  the  calculated  and  observed  data.  Both  calculated  and 
observed  trajectories  show  similar  patterns.  At  full-scale  depth  of  12.2  m  (i.e.,  40  ft),  the 
bomb  has  rotated  180°  such  that  the  bomb  is  moving  tail  first.  When  the  nose  reaches  a 
depth  of  12.2  m,  the  nose  horizontal  position  is  between  5.5  m  to  7.2  m.  Launch-13  had  a 
deviation  of  4.3°  from  the  vertical  water  entry  angle,  which  was  the  largest  deviation  in 
all  of  the  launches.  However,  comparison  of  the  trajectories  in  Launch-13  and  Launch-14, 
in  which  a  90°  water-entry  angle  was  obtained,  indicates  that  there  is  no  correlation 
between  impact  angle  and  trajectory  path.  Thus,  the  trajectory  motion  is  dominated  by 
instability  of  the  bomb  within  the  cavitated  region.  Different  from  the  Mk-84  bomb  with 
a  tail  section  and  no  fin  (or  two  fins),  the  Type-IV  bombs  are  never  moving  up  toward  the 
surface  (Fig.  18). 

Figs.  19-21  show  the  comparison  between  predicted  and  observed  time  evolutions 
of  the  horizontal  deviation  (x)  and  depth  position  (z)  of  Om,  bomb  speed  {U),  and  angles 
{a,/3,y  ).  The  predicted  values  of  these  variables  are  consistent  to  the  corresponding 
observed  values.  This  confirms  the  validity  of  the  semi-empirical  formulas  (53)-(55)  for 
drag/lift  and  torque  coefficients  (Q,  C/,  Cm).  The  three  launches  show  the  same 
interesting  results.  The  bomb  nose  reaches  the  12.2  m  depth  at  about  110  ms  after  water 
entry.  At  this  depth  the  bomb  nose  velocity  decreased  by  about  82%.  The  horizontal 
deviation  (x)  of  Om  increases  rapidly  at  first  and  then  slowly  with  time,  and  about  6  m 
from  the  entry  point  as  the  bomb  reached  the  depth  of  12.2  m. 
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8.  Tail  Section  Damage  Effects 

The  experiments  eondueted  with  different  tail  eonfigurations  were  performed  to 
determine  the  effeets  on  bomb  trajeetory  for  different  possible  postulated  damage  levels 
to  the  tail  seetion.  Beeause  the  tail  seetion  is  eomprised  mostly  of  internal  stiffeners  with 
an  external  skin,  it  may  be  weaker  than  the  warhead  seetion  and,  therefore,  may  be 
damaged  during  initial  water  entry  or  during  tail  slap  within  the  eavitated  region.  Fig.  22 
shows  the  maximum  measured  trajeetories  for  eaeh  tail  eonfiguration  for  a  nominal 
water-entry  veloeity  of  about  297  ms’',  i.e.,  Launoh-2  (Type-I),  Launoh-19  (Type-II), 
Launoh-17  (Type-III),  and  Launoh-13  (Type-IV). 

Table  2  summarizes  data  eomparison  of  the  overall  trajeetory  behavior  for  the 
different  tail  eonfigurations.  For  eaeh  trajeetory  parameter  we  show  the  value  assoeiated 
with  a  partieular  tail  eonfiguration  and  the  pereentage  differenee  eompared  with  a 
eomplete  bomb  having  a  tail  seetion  and  four  fins.  The  horizontal  nose  position 
signilioantly  inereases  with  inereased  levels  of  damage  to  the  tail  seetion.  For  a  model 
with  a  tail  seetion  and  two  fins,  no  fins,  and  no  tail  seetion,  the  horizontal  position  values 
inerease  by  120%,  259%,  and  575%,  respeetively.  With  regard  to  travel  time  at  12.2  m 
depth,  only  the  no-tail  eonfiguration  shows  a  signifieant  inerease  of  179%. 

9.  Conclusions 

(1)  A  new  dynamie-photographie  method  has  been  developed  to  determine  the 
drag/lift  and  torque  eoeffieients  (C^,  C/,  Cm)  of  a  fast-moving  rigid  body  in  the  water 
eolumn.  This  method  eontains  two  parts:  (1)  establishment  of  the  diagnostie  relationship 
between  (C*  C/,  Cm)  and  the  rigid-body’s  trajeetory  and  orientation,  and  (2)  data 
oolleetion  of  trajeetory  and  orientation  of  a  fast  moving  rigid  body  using  multiple  high- 


23 


speed  video  eameras  (10,000  Hz).  Using  the  digital  photographie  data,  semi-empirieal 
formulas  of  (Cj,  C/,  Cm)  versus  Reynolds  number,  attaek  angle,  and  rotation  rate  ean  be 
established.  The  eost  of  this  method  is  mueh  lower  than  the  traditional  method  using  the 
wind  tunnel  to  determine  (C*  C/,  Cm). 

(2)  To  demonstrate  the  feasibility  and  powerfulness  of  this  method,  an  experiment 
was  eondueted  with  1/12**'  seale  model  of  the  general  purpose  bomb  (Mk-84)  as  the  fast 
moving  rigid  body  at  a  6  m  deep  by  9  m  diameter  pool,  loeated  at  the  SRI’s  Corral 
Hollow  Experiment  Site.  A  gas  gun  was  used  to  shoot  the  1/12**'  model  Mk-84  bomb  into 
the  water  tank  with  veloeities  up  to  304  ms‘\  Four  types  of  Mk-84  model  bombs  were 
used  with  total  16  launehes  for  the  experiment:  warhead  with  tail  seetion  and  four  fins 
(Type-1),  with  tail  seetion  and  two  fins  (Type-H),  with  tail  seetion  and  no  fm  (Type-HI), 
and  with  no  tail  seetion  (Type-IV).  Among  them,  data  from  3  launehes  in  Type-lV  were 
used  to  get  the  semi-empirieal  formulas  for  (C^,  C/,  Cm).  The  rest  of  data  were  used  for 
verifieation. 

(3)  The  momentum  equations  and  moment  of  momentum  equations  were 
integrated  with  the  same  parameters  (sueh  as  the  density  ratio,  length,  radius,  the  eenter 
of  volume,  and  the  eenter  of  mass)  and  the  drop  initial  eonditions  (speed  and  orientation) 
as  in  the  observations  after  using  the  semi-empirieal  formulas  for  (C*  C/,  Cm).  Consistent 
between  ealeulated  and  observed  bomb  trajeetories,  orientations,  and  veloeities  show  the 
powerfulness  of  this  method. 

(4)  Both  ealeulated  (solving  dynamie  equations  with  the  semi-empirieal  formulas) 
and  experimental  data  show  similar  results.  The  eavitated  eolumn  generated  by  the  Mk84 
bomb  motion  eonsists  of  a  tapered  eone  that  has  a  maximum  full  seale  diameter  at  the  end 
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of  the  tail  section  of  about  0.9  m.  This  is  about  a  factor  of  two  larger  than  the  maximum 
bomb  diameter  of  0.45  m.  This  cavity  shape  was  about  the  same  for  all  of  the  initial 
water-entry  velocities  and  four  types  of  model  bombs. 

(5)  The  horizontal  nose  position  significantly  increases  with  increased  levels  of 
damage  to  the  tail  section.  For  a  model  with  a  tail  section  and  two  fins,  no  fins,  and  no 
tail  section,  the  horizontal  position  values  increase  by  120%,  259%,  and  575%, 
respectively.  With  regard  to  travel  time  at  12.2  m  (i.e.,  40-ft)  depth,  only  the  no-tail 
configuration  shows  a  significant  increase  of  179%. 

(6)  For  bomb  with  a  tail  section  and  four  fins,  its  trajectories  are  quite  stable 
without  oscillation  and  tumbling  no  matter  the  water  entry  velocity  is  high  or  low. 
Removal  of  two  fins  causes  the  bomb  to  eventually  make  180°  turn  and  travel  toward  the 
surface.  Although  having  a  similar  trajectory  pattern  (i.e.,  makingl80°  turn  and  traveling 
toward  the  surface),  the  removal  of  four  fins  shows  a  decrease  in  overall  trajectory 
stability  compared  with  a  model  with  a  tail  section  and  two  fins.  This  is  evident  primarily 
through  larger  horizontal  positions.  For  Mk-84  bomb  without  a  tail  section,  the  bomb  has 
rotated  180°  at  full-scale  depth  of  12.2  m  such  that  the  bomb  is  moving  tail  first.  Different 
from  the  Mk-84  bomb  with  a  tail  section  and  no  fin  (or  two  fins),  the  Mk-84  bombs 
without  a  tail  section  are  never  moving  up  toward  the  surface. 
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Table  1.  Summary  of  Mk-84  underwater  trajectory  experimental  matrix. 


Launch 

Number 

Model  Type 

Water-Entry 
Velocity  (m/s) 

Water- 

Entry 

Impact 

Angle  C) 

1 

I  (tail  with  4  fins) 

132 

89.2 

2 

I  (tail  with  4  fins) 

297 

90.0 

3 

I  (tail  with  4  fins) 

295 

88.8 

4 

I  (tail  with  4  fins) 

302 

88.5 

5 

I  (tail  with  4  fins) 

227 

88.0 

6 

I  (tail  with  4  fins) 

219 

89.0 

7 

I  (tail  with  4  fins) 

119 

88.2 

8 

II  (tail  with  2  fins) 

impacted 

sabot 

stripper 

plate 

9 

II  (tail  with  2  fins) 

impacted 

sabot 

stripper 

plate 

10 

II  (tail  with  2  fins) 

295 

90.0 

11 

II  (tail  with  2  fins) 

290 

90.0 

12 

II  (tail  with  2  fins) 

impacted 

sabot 

stripper 

plate 

13 

IV  (no  tail) 

297 

85.7 

14 

IV  (no  tail) 

301 

90.0 

15 

IV  (no  tail) 

301 

88.7 

16 

III  (tail  with  no  fin) 

304 

90.0 

17 

III  (tail  with  no  fin) 

297 

87.0 

18 

III  (tail  with  no  fin) 

291 

88.1 

19 

II  (tail  with  2  fins) 

297 

90.0 
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Table  2,  Horizontal  position  shift  and  travel  time  at  depth  of  12,2  m  (i,e,,  40  ft)  for 
Mk-84  warhead  with  different  tail  configurations. 


Model  Type 

Horizontal  Position 
Shift  (Unit;  m) 

Travel  Time 
(Unit:  ms) 

Tail  with 

4  Fins 

0.96 

46.2 

Tail  with 

2  Fins 

2.12 

46.2 

Tail  with 
Non  Fins 

3.46 

49.8 

No  Tail 

6.50 

129.5 
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Fig.  1 .  Position  vectors  r/,,  r^,  and  the  unit  vector  e. 
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Fig.  2.  Attack  angle  (a),  angles  center  of  volume  (ov),  center  of  mass  (om),  and 

drag  and  lift  forces  (exerted  on  Oy).  Note  that  a  is  distanee  between  Ov  and  Om  with 
positive  (negative)  value  when  the  direction  from  Oy  to  Om  is  the  same  (opposite)  as  the 
unit  veetor  e. 
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Fig.  3.  Illustration  of  unit  vectors  (e„,  e^,  e^,  e^^)with  ef  the  projection  of  e„  on  the 
x-y  plane. 
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Fig.  4.  Axial  and  cross-axial  velocity  (uf,  uf  ),  associated  hydrodynamie  forces  on  a 

pair  of  fins  (F/,  F/  )  and  the  distance  between  o/  and  Om  (i.e.,  cr^  )  with  positive 

(negative)  value  when  the  direetion  from  o/  to  Om  is  the  same  (opposite)  as  the  unit 
vector  e. 
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Fig.  5.  Illustration  of  H,  M^. 
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(a) 


(b) 


Fig.  6.  Photography  of  1/12*  scale  model  Mk-84  bomb;  (a)  warhead  with  tail  section 
and  four  fins,  and  (b)  sabot. 
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Fig.  7.  Overview  experimental  arrangement. 
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Fig.  8.  Two  HSV  images  for  Launch-3  (Type-I)  at  water  entry  velocity  of  295  ms"'; 
initial  water  entry,  (b)  t  =  22.8  ms,  and  (c)  t  =  44.4  ms. 


37 


(a)  -5/s 


Fig.  9.  Dependence  of  Cd  on  the  Reynolds  number  (Re)  and  attack  angle  ( «  )  for  three 
different  values  of  Q  ;  (a)  -5  s'\  (b)  0,  and  (c)  5  s'\ 
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(a)  £2=  -5/s 


Fig.  10.  Dependence  of  C/  on  the  Reynolds  number  (Re)  and  attack  angle  ( « )  for  three 
different  values  of  Q  ;  (a)  -5  s'\  (b)  0,  and  (c)  5  s'\ 
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0.15 


Fig.  1 1 .  Dependence  of  Cm  on  the  Reynolds  number  (Re)  and  attack  angle  {a). 
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(a)  (b)  (c)  (d)  (e)  (f)  (g) 


X(m)  X(m)  X(m)  X(m)  X(m)  X(m)  X(m) 


Fig.  12.  Comparison  between  predicted  and  observed  trajectories  for  Mk84  warhead  with 
tail  section  and  four  fins  (Type-1)  with  initial  water  entry  speed:  (a)  132  ms  (b)  297 
ms'\  (c)  295  ms'^  (d)  302  ms'\  (e)  227  ms'^  (f)  219  ms'^  and  (g)  119  ms'\ 
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Fig.  13.  Two  HSV  images  for  Launch- 1 1  (Type-II)  at  water  entry  velocity  of  290  ms'*; 
(a)  initial  water  entry,  (b)  t  =  21.6  ms,  (c)  t  =  48.0  ms,  (d)  t  =  75.6  ms,  (e)  t=  1 16.4  ms, 
and  (f)  t  =  344.4  ms. 


42 
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Y(m)  Y(m)  Y(m) 

Fig.  14.  Comparison  between  predieted  and  observed  trajeetories  for  Mk84  warhead  with 
tail  section  and  two  fins  (Type-11)  with  initial  water  entry  speed;  (a)  295  ms  (b)  290 
ms'\  and  (c)  297  ms’'. 


(a)  (b)  (c)  (d) 


Fig.  15.  Two  FISV  images  for  Launch-17  (Type-Ill)  at  water  entry  velocity  of  298  ms"': 
(a)  initial  water  entry,  (b)  t  =  22.8  ms,  (c)  t  =  55.2  ms,  (d)  t  =  99.0  ms,  (e)  t  =  21 1.2  ms, 
and  (f)  t  =  376.2  ms.  Note  that  for  time  longer  than  99.0  ms,  only  one  FISV  camera  got 
the  pictures. 
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(a)  Initial  U  =  304m/s  (b)  Initial  U  =  299m/s  (c)  Initial  U  =  291  m/s 


Fig.  16.  Comparison  between  predicted  and  observed  trajectories  for  Mk84  warhead  with 
tail  section  and  no  fin  (Type-111)  with  initial  water  entry  speed;  (a)  304  ms'^  (b)  298  ms‘\ 
and  (c)  291  ms’'. 
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Fig.  17.  Two  FISV  images  for  Launch-13  (Type-IV)  at  water  entry  velocity  of  296  ms"'; 
(a)  initial  water  entry,  (b)  t  =  30.0  ms,  (c)  t  =  51.6  ms,  (d)  t  =  155.4  ms,  and  (e)  t  =  418.2 
ms. 
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(a)  Initial  U  =  296m/s 


Experiment 


3D  Modei 


Fig.  18.  Comparison  between  predicted  and  observed  trajectories  for  Mk84  warhead  with 
no  tail  section  (Type-lV)  with  initial  water  entry  speed;  (a)  296  ms  (b)  301  ms'\  and 
(c)  301  ms'\ 
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U  (m/s)  Z  (m) 


Fig.  19.  Time-evolutions  between  predieted  (solid)  and  observed  (dotted)  for  Launoh-13 
(a)  horizontal  position  (x)  of  Om,  (b)  depth  position  (z)  of  Om,  (c)  bomb  speed  (U),  (d) 
angle  y ,  (e)  angle  ,  and  (f)  attaek  angle  a  . 


U  (m/s)  Z  (m) 


(c)  (f) 


Fig.  20.  Time-evolutions  between  predieted  (solid)  and  observed  (dotted)  for  Launoh-14: 
(a)  horizontal  position  (x)  of  Om,  (b)  depth  position  (z)  of  Om,  (c)  bomb  speed  (U),  (d) 
angle  y ,  (e)  angle  P ,  and  (f)  attaek  angle  a  . 


Time  (s)  Time  (s) 


Fig.  21.  Time-evolutions  between  predieted  (solid)  and  observed  (dotted)  for  Launoh-15; 
(a)  horizontal  position  (x)  of  Om,  (b)  depth  position  (z)  of  Om,  (c)  bomb  speed  (U),  (d) 
angle  y ,  (e)  angle  P ,  and  (f)  attaek  angle  a  . 
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